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ABSTRACT 

There is now strong evidence that some stars have been born with He mass fractions 
as high as y « 0.40 (e.g., in uj Centauri). However, the advanced evolution, chemical 
yields, and final fates of He-rich stars are largely unexplored. We investigate the 
consequences of He-enhancement on the evolution and nucleosynthesis of intermediate- 
mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 Mq with a metallicity 
of Z = 0.0006 ([Fe/H] « —1.4). We compare models with He-enhanced compositions 
(Y = 0.30,0.35,0.40) to those with primordial He {Y = 0.24). We find that the 
minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores 
decreases from above our highest mass of 6 Mq to ~ 4-5 Mq with Y = 0.40. We 
also model the production of trans-Fe elements via the slow neutron-capture process 
(s-process). He-enhancement substantially reduces the third dredge-up efficiency and 
the stellar yields of s-process elements (e.g., 90% less Ba for 6 Mq, Y = 0.40). An 
exception occurs for 3 Mq , where the near-doubling in the number of thermal pulses 
with Y = 0.40 leads to ~ 50% higher yields of Ba-peak elements and Pb if the ^^G 
neutron source is included. However, the thinner intershell and increased temperatures 
at the base of the convective envelope with Y = 0.40 probably inhibit the neutron 
source at this mass. Future chemical evolution models with our yields might explain 
the evolution of s-process elements among He-rich stars in w Gentauri. 

Key words: stars: AGB and post-AGB stars — stars: evolution — nuclear reactions, 
nucleosynthesis, abundances 


1 INTRODUCTION 


Stars that evolve through the asymptotic giant branch (AGB) 
play a crucial role in the chemical evolution of stellar popula¬ 
tions and galaxies. These stars experience a complex sequence 
of He-shell instabilities, mixing, and mass loss processes that 
combine to eject material that has undergone H and He 
burning, and enrichment in heavy elements produced by the 
slow neutron-capture process (s-process). Understanding the 
chemical contribution by stars to the interstellar medium is 
a prerequisite for building models of chemical evolution (e.g., 
Chiappini et al.|2001 Kobayashi et al.||201l K 

The most important parameter governing stellar evo¬ 
lution is the initial mass, with a secondary role played by 
chemical composition. Chemical composition typically refers 
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to the mass fraction of metals (Z), but the mass fraction 
of He (y) also has major consequences for stellar evolution. 
However, the stellar evolution and nucleosynthesis that take 
place during the AGB phase are relatively unexplored for 
He-rich initial compositions. For a detailed introduction to 
AGB evolution and nucleosynthesis with normal He content, 
we refer the reader to the reviews by Herwig (20051 and 


Karakas & Lattanzio (20141. 


There is strong evidence that some stars have been born 
with substantial He enrichments above the primordial He 
mass fraction of F « 0.24 predicted from big bang nucleosyn¬ 
thesis. The most direct method to measure He abundances is 
to use spectroscopy, however the application to He is severely 
hindered by the lack of He lines in cool stars, and the effects 
of gravitational settling in hot stars. Still, direct detections of 
He have been made, such as those by |Pasquini et aL] ( |2011[ ), 
who used the He I 10830 A line in two stars with different 
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Na and O abundances in NGC 2808 to derive a difference of 
AY ^ 0.17. The finding of He-enhanced second-generation 


stars in NGC 2808 is also supported by Marino et al. (20141 
with a sample of 96 horizontal-branch stars, although they 
report a smaller average enhancement of AY = 0.09 ± 0.06. 


(hereafter VD09), which have Z = 0.0006 and [a?/Fe] = 0.4. 
The models of VD09 are calculated with the Full Spectrum 
of Turbulence (FST, Canuto fc Mazzitelli||199r I treatment of 


The same He 1 10830 A line was also used by [Pupree fc 


convective mixing and energy transport rather than the more 
common Mixing Length Theory (MLT, Biermann||1948 1. 


Avrett (20131 to infer a variation of AY ^ 0.17 between red 


giants in uj Centauri. 

A second, less-direct line of evidence for He-rich stars 
is the requirement of high He abundances to reconcile theo¬ 
retical isochrones with the photometry of clusters in colour- 
magnitude diagrams (CMDs). For example, the split main 
sequence of uj Centauri can be explained by multiple pop- 


The stellar yields of He-rich models, including those pre¬ 
sented in this study will build a foundation for future work 
in understanding the chemical evolution of He-rich environ¬ 
ments. Our comparison with other models in the literature 
allows us to estimate how sensitive our predictions for He- 
rich stellar evolution and yields are to uncertain modelling 
assumptions. 


ulations that vary in both metallicity and He abundance 

and stellar modelling assumptions. In Section [^ we present 

Bedin et al.|2004 Norris|2004 Piotto et al.|2005 1. King et al. 

our stellar evolutionary models. In Section [^we present 

2012P report a best fit to the blue main sequence with a He 

the nucleosynthesis models and explore the impact of He- 


abundance of T = 0.39 ± 0.02. The high He abundance also 
persists when additional features of the CMD are incorpo¬ 
rated in a simultaneous fit, such as the work of|Joo fc Lee| 


(20131, who infer Y of up to 0.39 ± 0.02 for oj Centauri and 


up to 0.32 ± 0.04 for M22. 

With relatively few theoretical studies of the detailed 
AGB evolution and nucleosynthesis of He-rich models (|Camp-| 


bell et al.|2013| [Charbonnel et al.|2013| [Karakas et ai.||2014| 

Karakas|2014 1, the effects of high He abundances on chem¬ 

ical evolution are not well understood. An important set 
of observational clues are the s-process abundances in uj 
Centauri. In this cluster, [Ba/FeQ increases with metallicity 
up to about [Fe/H] = —1.4, above which [Ba/Fe] remains 
roughly constant ( [Norris fc Da Costa||199^ ). Interestingly, 
this plateau occurs near the metallicity of the most He-rich 
population ( Joo fc Lee|2013 1, and one potential explanation 
is that Ba yields are lower in He-rich stars. Chemical evolu¬ 
tion might play out similarly in other clusters such as M22, 


which follows the same rise of [Ba/Fe] (Da Costa & Marino 


20111 but where a truncation of the chemical evolution at a 


metallicity below [Fe/H]= —1.4 prevents the possibility of a 
similar s-process plateau. 

Karakas, Marino &: Nat^ ( [2014 ) (hereafter KMN14) 


present stellar yields of He-enhanced 1.7 Mq and 2.36 Mq 
AGB models with a metallicity of [Fe/H] « —1.4. This metal¬ 
licity is appropriate for the most metal-poor of the He-rich 
populations in lo Centauri ( Joo fc Leej 20131. Among the 
findings of KMN14 is a reduction in the third dredge-up 
mass with increasing initial He abundance. Largely as a re¬ 
sult of the less efficient dredge up, their Ba yields decrease 
by roughly 50 per cent with Y = 0.40 in comparison with 
primordial-He models at the same initial mass. 

In this work, we extend the study of He-enhanced AGB 
nucleosynthesis by KMN14 with new intermediate-mass mod¬ 
els of 3, 4, 5, and 6 Mq at the same metallicity of Z — 0.0006 
and with similar input physics. For our primordial He models, 
we make a comparison with the models of jStraniero, Cristalloj 
jfc Piersanti| ( |201^ (hereafter SCP14), which have Y = 0.245, 
Z = 0.0003 (before a-enhancement, [Fe/H] = —1.7), and 
[a/Fe] = 0.5. We also compare our evolutionary models with 
the primordial-He models of [Ventura fc D’Anton^ (|2009[| 


^ We use the standard spectroscopic notation, [A/B] = 

logioC-^A/Ns) - logiQ(A®/A''®), where Na and Nb are abun¬ 
dances by number and © denotes the solar abundance. 


our results in the context of other studies in the literature 
and consider the implications for He-rich chemical evolution. 


2 COMPUTATIONAL METHOD 

We present evolutionary sequences for initial masses of 3, 4, 
5, and 6 Mq with a global metallicity of Z = 0.0006 ([Fe/H] 
~ —1.4 and [a/Fe]= 0.0), and initial He mass fractions of 
Y = 0.24,0.30,0.35, and 0.40. 

We use the same input physics and stellar evolutionary 
code as KMN14 (the Mount Stromlo Stellar Structure 
Program), except that we include the opacity treatment 
of [Fishlock et ^ (20141, which accounts for changes to 
the C abundance and the C/O ratio at the surface. This 
opacity treatment improves the accuracy for models that 


experience hot bottom burning (discussed in Section 3.41. 
At low-temperatures (T < 10,000 K), we use the molecular 
opacity tables from AESOPUS (Marigo & Aringer 20091. 


At higher temperatures, we use OPAL tables ( Iglesias fc 
Rogers[[l996 1 that have been updated to a Lodders (20031 
scaled-solar composition. 

The initial abundances of the models are scaled solar 
( Asplund et al.|2009 1 with a global metallicity of Z = 0.0006. 
With Z fixed and a chosen value of Y, the H mass fraction 
(V) is determined byX = l — V - Z. A consequence of 
keeping the mass fraction of metals (including Fe) constant 
and exchanging mass between H and He is that the initial 
ratio of Fe to H varies slightly between our models of different 
He abundance. Specifically, the initial [Fe/H] is —1.41 in the 
models for Y = 0.24, and —1.31 in the models for Y = 0.40. 

The search for stable convective borders uses the method 
of Lattanzio (19861 (see also Frost fc Lattanzio|1996 1, and we 
assume no other form of extra mixing at convective bound¬ 
aries. The convective velocities are calculated using MLT 
with a mixing length parameter of a = 1.86. Although the 
constancy of a is a standard assumption, empirical and theo¬ 
retical studies have found a variation of this parameter with 
stellar evolution ( Lebzelter fc Wood[[2007 Magic & Asplund 
20141. Larger values of a can increase the depth of the third 


dredge-up (e.g., Boothroyd & Sackmann 19881, which in 
turn typically increases the yields of C-l-N-l-0 and s-process 
elements ( Cristallo et al.|2009 20111. 

We do not include mass loss on the first red giant branch 
(RGB). This assumption is reasonable (in the context of 
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stellar yield predictions) given recent results using Kepler 
data, which indicate that the total mass loss during the RGB 


Hollowell &: Iben|1989K Further, the important effect of ro- 


tationally induced mixing on the neutron release in the '^C 


is relatively small (AM = 0.09 ±0.03 ±0.04 Mq, Miglio et al. 


20121, at least near solar-metallicity. Many of our models 


pocket has been investigated by Herwig et al. (20031, Siess 
et al. (20041, and Piersanti et al. (20131. 


do not reach the RGB and for the models that do, the 
inclusion of mass loss of the this order would have a minimal 
impact on our stellar yield predictions, which are subject 
to much larger uncertainties related to mass loss during the 


AGB phase (Stancliffe & Jeffery 20071. For mass loss on 
the AGB, we use the prescription given by Equations 1 to 4 


of Vassiliadis & Wood (1993), which includes a superwind 


phase beginning at pulsation periods above approximately 
500 days. 

The output of the stellar evolution code (including the 
temperatures, densities, convective boundaries, and convec¬ 
tive velocities) is used as input to a post-process code that 
performs the detailed nucleosynthesis calculations with a 
large network of 320 species up to Po. The two-pass method 
used in this work is different to the models of SCP14 cal¬ 
culated with the FUNS code (a descendant of FRANEC, 
Chiefli fc Straniero|1989 1, which solves the stellar structure 


equations simultaneously with the full network of nuclear 
reactions. 

Another difference is that our post-processing code 
solves simultaneously the abundance changes wrought by 
mixing and burning. To mix the convective regions in 
our nucleosynthesis models, we use the two-stream, time- 


The nature of the transition away from pockets with 
increasing initial mass is also highly uncertain, although it 
is very likely connected with the shrinking of the He-rich 
intershell and the higher temperatures at the base of the 
convective envelope with larger core masses. 

For example, the models of SCP14 assume that the con¬ 
vective velocity declines exponentially with distance beyond 
the formal convective boundary, with a free parameter de¬ 
scribing the length scale of the velocity decline. This method 
is applied only to convective boundaries for which the veloc¬ 
ity is discontinuous, which occurs at the base of the envelope 
during third dredge-up episodes ( Cristallo et al.|2009 ). The 
convective velocity is incorporate d into a non-diffusive mix¬ 
ing scheme ( Straniero et al.|2006 l and results in pockets 
that smoothly decrease in size (and neutron production) with 
increasing initial mass. 


Goriely & Siess (20041 show that a diffusive treatment of 


convective-boundary mixing can reduce the amount of neu¬ 
tron production by pockets (due to enhanced production 
dependent method devised by |Cannoii| (|1993|). An upward of ^"‘N) when the third dredge-up takes place with temper¬ 
atures above 40 MK, and can totally inhibit the s-process 
with temperatures above 70 MK. As shown in Table the 
maximum temperatures at the base of the envelope during 
third dredge-up for our 3 Mq models with Y between 0.24 
and 0.35 occupy this transition range of temperatures. For 
the 3 Mq model with Y — 0.40 and the higher-mass models 
with all He abundances, temperatures above 70 MK during 
third dredge-up would likely prevent a significant s-process 
production by pockets. For this reason, a PMZ leading to 
the formation of a pocket is not included in our models 
with M ^ 4 Mq. For the initial mass of 3 Mq, we construct 
nucleosynthesis models with and without a PMZ at each He 
abundance. 


and a downward-moving stream are treated separately, with 
the mass-flow rate at each mass shell calculated from the 
convective velocities in the evolutionary model. It is assumed 
that the material from an upward moving cell below flows 
into the upward moving cell above, and likewise the material 
in a downward moving cell flows into the downward moving 
cell below, as well as horizontally mixing from one stream to 
the other. 


2.1 The inclusion of pockets 


For stars with initial masses < 3 Mq, the dominant neutron 
source dr iving s-process nu cleosynthesis is the ^^C(a,n)^®0 
reaction ( Bu^o et al.|l99'9 l, while the ^^Ne(Q:, n)^®Mg reac¬ 
tion becomes more important at higher initial masses ( |Garcfa-| 
Hernandez et al.|2013 1 . For our 3 Mq models which are near 


the mass of this transition, the treatment of ^^G-pocket for¬ 
mation is crucial for making accurate predictions of the stellar 
yields and surface abundances of neutron-capture elements. 

In current theoretical models, the ^^C(a,n)^®0 reaction 
occurs in a thin layer of that results from the partial 
mixing of protons from the base of the envelope into the 


In the absence of a deep physical understanding of the 
PMZ, we apply the simple parameterised treatment used by 
Lugaro et al. (20041 (similar to Goriely fc Mowlavi|2000| in 


C-rich intershell region during third dredge-up (Straniero 


our nucleosynthesis post-processing models. For models that 
include a PMZ, we insert an exponential (in mass coordinate) 
profile of protons below the envelope at the deepest extent of 
each third dredge-up episode. The proton abundance starts 
at the envelope value and decreases by a factor of 10^ through 
a depth given by the parameter Mpmz. The value of Mpmz 
is uncertain, although a variety of constraints have been 


in N, which efficiently absorb free neutrons. The physical 
mechanism that produces the partially mixed zone (PMZ) 
has not been conclusively identified, but plausible candidates 
include convective-boundary mixing ( Herwig|2000 Cristallo 


tured by where they activate the CN-cycle reactions 

derived from observations of AGB stars (Abia et al.||2002|. 

^^C(p, 7 )^®N(/?+i/)^®C. The number of protons mixed into 

C-enhanced metal poor stars (Izzard et al.||2009 

Bisterzo 

the region must be low because further proton capture results 

et al.|2012| |Lugaro et al.|2012||, planetary nebulae 

Shingles 


(Bonacic Marinovic et al. 2007 De Smedt et al. 2012 I. KMN14 


et al.|20^), rotational mixing (Langer et al.||19^ Herwig 


fc Langer||2001[), gravity-wave driven mixing ([Denissenkov 


& Tout 2003), and semiconvection (Iben & Renzini 1982 


presented 1.7 Mq and 2.36 Mq models with a range of PMZ 
sizes, but for comparison we only use their results for a 
PMZ of 10~^ Except where stated otherwise, our 3 

Mq nucleosynthesis models include a PMZ with Mpmz = 
10"^ Mq. 
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Table 1. Evolutionary lifetimes, first and second dredge-up masses, and final model properties. We include the initial stellar mass 
(Mini), fhe initial He mass fraction (E), the core H-burning lifetime (rms), the core He-burning lifetime (rcoreHe)? the AGB lifetime 
measured from the exhaustion of core He until the end of calculations plus the estimated time to lose the remaining envelope (ragb), 
the total stellar lifetime (Tsteiiar), the innermost mass layer reached by first and second dredge-up (Mfdu the final 

stellar mass the final core mass the final envelope mass the final mass-loss rate and the class of 

white dwarf remnant formed. Final refers the state at the end of our calculations, which cease due to convergence problems before 
the envelope has been fully ejected. 


Mini 

Y 

rms 

^coreHe 

■Tagb 

"^stellar 

MpDU 

MsdU 

Mflnal 

Mflnal 

Mfl-> 

Mflnal 

Remnant 



[Myr] 

[Myr] 

[Myr] 

[Myr] 

[Mq] 

[Mq] 

[M©] 

[Mq] 

[Mq] 

[10-5 Mq] 



0.24 

213.2 

59.6 

0.81 

290.1 

1.79 

0.786 

1.462 

0.801 

0.661 

2.64 

CO WD 

3 Mq 

0.30 

161.2 

45.5 

0.29 

219.5 

2.28 

0.844 

1.182 

0.858 

0.324 

3.06 

CO WD 

0.35 

126.1 

35.9 

0.23 

172.5 

2.94 

0.863 

1.628 

0.880 

0.748 

3.39 

CO WD 


0.40 

97.5 

29.0 

0.20 

134.8 

- 

0.888 

1.641 

0.914 

0.727 

3.98 

CO WD 


0.24 

112.4 

26.6 

0.97 

146.3 

- 

0.862 

1.513 

0.886 

0.626 

3.56 

CO WD 

4 Mq 

0.30 

85.6 

21.4 

0.61 

112.4 

- 

0.887 

1.462 

0.910 

0.552 

4.23 

CO WD 

0.35 

67.8 

17.5 

0.36 

89.5 

- 

0.913 

1.625 

0.938 

0.687 

4.12 

CO WD 


0.40 

53.2 

14.5 

0.20 

71.0 

- 

0.951 

1.760 

0.984 

0.776 

4.73 

CO WD 


0.24 

71.1 

14.6 

0.67 

89.3 

- 

0.913 

1.671 

0.940 

0.731 

4.02 

CO WD 

5 Mq 

0.30 

54.8 

11.9 

0.33 

69.4 

- 

0.951 

1.597 

0.976 

0.621 

4.32 

CO WD 

0.35 

43.7 

10.3 

0.18 

56.0 

- 

0.998 

1.603 

1.025 

0.578 

4.92 

CO WD 


0.40 

34.5 

8.9 

0.11 

45.0 

- 

1.073 

1.911 

1.125 

0.786 

5.99 

CO(Ne) WD 


0.24 

50.0 

9.5 

0.39 

61.4 

- 

0.996 

1.569 

1.021 

0.549 

4.43 

CO WD 

6 Mq 

0.30 

38.8 

8.0 

0.13 

48.0 

- 

1.058 

1.624 

1.080 

0.544 

5.09 

CO(Ne) WD 

0.35 

31.1 

6.7 

0.07 

38.8 

- 

1.140 

1.632 

1.168 

0.464 

6.04 

ONe WD 


0.40 

24.9 

6.0 

0.04 

31.6 

- 

1.234 

2.157 

1.258 

0.900 

9.15 

ONe WD 


Table 2. Structural properties relevant to the thermally pulsing AGB phase. We include the initial stellar mass (Mini), the initial He mass 
fraction (y), the H-exhausted core mass at the first thermal pulse (Mj^^), the number of calculated thermal pulses (TPs), the dredge-up 
efficiency A averaged over all TPs ((A)), the maximum dredge-up efficiency (Amax), the total TDU mass (Mi^^^), the maximum extent of 
the pulse-driven convective zone (M™^^), the maximum temperature in the He-burning shell the maximum temperature in 

the H-burning shell the maximum temperature at the base of the convective envelope during dredge-up or during interpulse 

H-burning and the average interpulse period {(rip)), and the total duration for which T^ce > 50MK (thbb)- 


Mini 

Y 

mTpi 

[M©] 

TPs 

(A> 

Amax 

[M©] 

Mmax 

pdcz 

[10-3 M©] 

rpmax 

He —shell 

[MK] 

'X'max 
-^H-shell 

[MK] 

rpmaXjdup 

-*-bce 

[MK] 

rpinaXjip 

^ bee 

[MK] 

(Pip) 

[Kyr] 

ihbb 

[Kyr] 


0.24 

0.786 

26 

0.94 

1.00 

0.160 

9.39 

361 

82 

57 

43 

27.71 

0 

3 Mq 

0.30 

0.847 

23 

0.86 

0.97 

0.070 

5.70 

338 

88 

63 

54 

11.61 

82 

0.35 

0.867 

26 

0.79 

0.94 

0.054 

4.36 

340 

92 

64 

63 

6.87 

86 


0.40 

0.892 

48 

0.70 

0.91 

0.065 

3.00 

344 

98 

78 

86 

3.63 

120 


0.24 

0.864 

82 

0.92 

0.99 

0.281 

4.78 

350 

90 

74 

84 

9.91 

683 

4 Mq 

0.30 

0.890 

82 

0.88 

1.00 

0.186 

3.48 

354 

95 

79 

88 

5.60 

382 

0.35 

0.918 

75 

0.81 

0.93 

0.114 

2.59 

355 

100 

83 

93 

3.41 

216 


0.40 

0.956 

84 

0.62 

0.87 

0.058 

1.70 

349 

109 

90 

101 

1.62 

122 


0.24 

0.915 

121 

0.90 

0.97 

0.254 

3.02 

360 

99 

82 

95 

5.08 

558 

5 Mq 

0.30 

0.955 

104 

0.84 

0.98 

0.137 

2.13 

361 

105 

89 

99 

2.85 

277 

0.35 

1.002 

106 

0.70 

0.88 

0.069 

1.36 

368 

115 

97 

109 

1.38 

139 


0.40 

1.077 

261 

0.26 

0.64 

0.018 

0.48 

356 

132 

118 

126 

0.34 

91 


0.24 

0.998 

136 

0.81 

1.01 

0.140 

1.82 

384 

111 

96 

108 

2.39 

315 

6 Mq 

0.30 

1.061 

130 

0.71 

0.95 

0.058 

0.90 

381 

123 

108 

118 

0.94 

123 

0.35 

1.143 

227 

0.41 

0.73 

0.021 

0.37 

382 

136 

124 

131 

0.27 

65 


0.40 

1.236 

363 

0.20 

0.30 

0.005 

0.10 

353 

134 

122 

129 

0.08 

32 


3 STELLAR EVOLUTION MODELS 

The evolutionary properties of our models and their depen¬ 
dencies on the initial He abundance are reported in Tables 
and and discussed below. 


3.1 Pre-AGB evolution and stellar lifetimes 

Figure [T| shows the evolutionary tracks of the 3 Mq and 6 
Mq models with Y = 0.24 and 0.40 in a Hertzsprung-Russell 
diagram. In He-enhanced models, the initial compositions 














He-rich AGB evolution and nucleosynthesis 5 



Figure 1. Evolutionary Hertzsprung-Russell tracks of 3 and 6 
Mq models at Y = 0.24 (black lines) and Y = 0.40 (orange lines) 
from the main sequence to the beginning of the thermally pulsing 
AGB phase. Line styles indicate the nuclear burning stages: core 
H burning (solid), shell H burning (dot dashed), core He burning 
(dashed), and shell He and H burning (dotted). 


have a higher mean molecular weight (as H is exchanged 
for He), which leads to higher luminosities during the core 
H-burning phase ( Roeser|1975 1. Compared with models for 
Y = 0.24, the models for Y = 0.40 are roughly twice as 
luminous during the main sequence. They are also more 
luminous at each corresponding stage of evolution due to the 
larger core masses of He-enhanced models. 

As shown by the H-burning time-scales and total stellar 
lifetimes in Table the increased H-burning rate and the 
reduced H fuel available in He-enhanced models lead to 
much shorter core H-burning lifetimes and shorter stellar 
lifetimes overall. This has important implications for chemical 
evolution, as He-enhanced stars will process and eject their 
enriched material much earlier than primordial-He stars at a 
given initial stellar mass. With Y = 0.40, the stellar lifetimes 
of the models are reduced by roughly 50 per cent compared 
to models with primordial He abundance. The shortening of 
stellar lifetimes by a similar factor with Y = 0.40 is also found 
by Charbonnel et al. ( 2013[ ) for low-mass models between 
0.66 Mq and 0.82 Mq at [Fe/H] = —1.56, and KMN14 for 
1.7 Mq and 2.36 Mq models with the same metallicity of 
Z — 0.0006 as the models presented here. 


3.2 The first and second dredge-up 

After the exhaustion of H in the core, the onset of H burning 
in a shell around the core causes the outer layers to expand 
and cool as a star begins to ascend the red giant branch. At 
this stage, the envelope convection zone can move inwards in 
mass and dredge up the products of partial H burning, which 
increases the surface abundances of and and de¬ 
creases the abundance of Dredge-up at this evolutionary 
stage is known as first dredge-up (FDU). 

Table shows the innermost mass layer reached by the 


convective envelope during FDU for models in which FDU 
occurs. The depth of FDU in the 3 Mq models decreases 
with increasing He abundance up to Y = 0.40, for which 
FDU does not take place. This is because these stars ignite 
He in their cores before their envelopes expand sufficiently to 
reach the first giant branch (Figure]^. We find no FDU for 
initial masses of 4 Mq and above at any He abundance. For 
the Y = 0.24 models, the lack of FDU with initial masses 
of 4 Mq and higher is consistent with the slightly higher- 
metallicity (Z = 0.001) models of Fishlock et al. (20141, for 
which FDU ceases between initial masses of 3.5 and 4.0 Mq. 

Following the core He-burning phase, the stars change 
to a dual shell-burning structure and begin their ascent of 
the AGB. As the He-shell is ignited and the H-shell becomes 
extinguished for the first time, the inner boundary of the 
convective envelope moves inwards in mass and dredges 
up the products of complete H burning, which increases 
the surface abundances of '^He and This is the second 
dredge-up (SDU). 

Table shows the innermost mass later reached by the 
SDU for each of our models. The SDU is deeper than the 
FDU in every model with FDU, and the depth of the SDU 
decreases with increasing He abundance. However, due to 
the increased core mass in the 6 Mq model with Y = 0.40, a 


“corrosive SDU” takes place (Gil-Pons et al.|2013 Doherty 


et al.|20T4(. In this model, the inner edge of the convective 


envelope reaches below the top of the CO core, which dredges 
up C and O to the surface. 


3.3 The thermally pulsing AGB and third 
dredge-up 

During the AGB phase, thermal pulse cycles are driven by 
thermal instabilities of the thin He-burning shell (|Iben|197^ 


Herwig 2005 The rapid flash burning of the He-shell releases 


an enormous amount of energy and causes a pulse-driven 
convective zone to form, which homogenises the composition 
of the He-intershell. As the energy moves upwards and ex¬ 
pands the star, the H-burning shell is effectively extinguished. 
With the H-burning barrier removed, the envelope convection 
zone can penetrate into the He-intershell in an event known 
as a third dredge-up (TDU) episode ( Iben|1975 Sackmann 
|1980b[ ). 

To show how the thermal-pulse cycle is altered by He- 
enhancement, we first measure the average time between 
thermal pulses (the interpulse). The left panel of Figure 
shows the average interpulse period versus the average 
core mass during the AGB for our models from 3 to 6 Mq 
with He mass fractions of Y = 0.24, 0.30, 0.35 and 0.40. The 
interpulse time-scale is mainly a function of the H-exhausted 
core mass | Paczyhsid]|1974 Christy-Sackmann & Paczyriski 


1979 


19751, which increases with the He content ([Becker fc Iben 


Lattanzio 19861. However, the average interpulse period 


of a He-enhanced model is lower than that of a primordial-He 
model with a higher initial mass and the same average core 
mass. The shorter interpulse time with He-enhancement for 
a given core mass is connected with the increased H-burning 
rate. With Y = 0.40, the average interpulse period is up 
to 50 per cent shorter than with Y = 0.24 and the same 
average core mass. The right panel of Figure shows that 
a similar, but marginally tighter relation exists when the 
average core mass is replaced by the average He-intershell 
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Figure 2. Left: The average interpulse period versus the average H-exhausted core mass during the thermally pulsing AGB phase. Right: 
The mean interpulse period versus the mean of the maximum extent of the pulse-driven convective zone (a proxy for the He-intershell 
size). Included on the plot are the 3, 4, 5, and 6 Mq models (from left to right) with He mass fractions of Y = 0.24 (black squares), 0.30 
(green diamonds), 0.35 (pink circles), and 0.40 (orange triangles). 


mass (approximated by the maximum extent of the pulse- 
driven convective zone). 

Figure shows the H-burning luminosity and the tem¬ 
perature and luminosity of the He-burning shell as a function 
of time from just prior to the first thermal pulse until after 
the second thermal pulse (of the Y = 0.24 case). To illustrate 
how a He-rich initial composition alters the behaviour of 
the burning shells and reduces the interpulse time-scale, the 
models shown are the 3 Mq model with Y = 0.24, the 3 Mq 
model with Y — 0.35, and the 4 Mq model with Y — 0.24. 
These three models enable us to separately compare the effect 
of increasing the He mass fraction and increasing the initial 
mass, which each individually lead to a larger core mass. 
The 3 Mq, Y = 0.35 model and the 4 Mq, Y = 0.24 model, 
which have very similar core masses (0.865 ± 0.002 Mq) and 
intershell masses (5 x 10“® ^q) the time of their hrst 
thermal pulses also exhibit a very similar time evolution of 
their He-shell temperature and luminosity. However, even 
with near-identical core masses, the onset of the next thermal 
pulse takes place sooner in the He-enhanced model, which 
has a higher H-burning luminosity. 

The results in Tables [T] and [2] are derived from our evo¬ 
lutionary sequences which end due to convergence difficulties 
before the H envelopes are completely ejected ( Sweigart|1999 


Lau et ^|2012 1. The AGB lifetimes have been calculated 


from the time spent on the AGB up to the end of our calcula¬ 
tions, plus an estimate of the time left to lose the remaining 
H envelope. This is estimated by dividing the envelope mass 
at the last computed model by the average mass-loss rate 
over the last few thousand models. For the number of thermal 
pulses neglected by our calculations due to the early termi¬ 
nation of the AGB, an upper limit is obtained by dividing 
the envelope ejection time by the average interpulse time. 


The number of additional thermal pulses is less than 15 per 
cent of the numbers given in Table for all cases except for 
the 6 Mq, Y = 0.40 model, which could experience up to 
~ 35 per cent more thermal pulses beyond the end of our 
calculations. 

The total number of thermal pulses experienced by each 
model depends on both the interpulse period and the time 
spent in the AGB phase, which is controlled by the mass- 
loss rate. The number of thermal pulses has implications 
for the stellar yields, as a smaller number of thermal pulses 
means that there will be fewer neutron-producing events and 
typically fewer TDU episodes. For models with hot bottom 
burning (|Iben|1975| |Renzini fc Voli||1981| [Boothroyd &: Sadc^ 


mann[|19^ Lattanzio||1992 1, the length of time spent in the 

AGB phase also has an impact on the surface abundances 
and yields of elements involved in proton-capture reactions 
at the base of the envelope. In the He-rich AGB models, the 
mass-loss rates are higher due to higher luminosities and gen¬ 
erally longer pulsation periods, and thus they have shorter 
AGB lifetimes. Our models with Y = 0.40 have shorter AGB 
lifetimes compared with Y = 0.24 by a factor of four at 3 
Mq, and by a factor of ten at 6 Mq, although these results 
are dependent on the mass loss prescription. 

The models presented here for Y = 0.24 experience 
significantly more thermal pulses than the models of SCP14 
with the same initial mass. For 3, 4, 5, and 6 Mq, we find 26, 
82, 121, and 136 thermal pulses, compared with 15, 23, 35, 
and 72 for SCP14, and 75, 91, 103, and 116 for VD09. The 
smaller number of thermal pulses in the models of SCP14 is 
largely due to their modification of the | Vassiliadis fc Wood| 
( 1993[) mass-loss prescription to include a minimum rate of 


10 MQ/yr at log P < 2.7 (described in 


Straniero et al. 


2006[). VD09 use the Bloecker (19951 mass-loss prescription 


























He-rich AGB evolution and nucleosynthesis 7 


3 Mg Y=0.24, Mj"”! =0.786 Mg 

3 Mg Y=0.35, =0.867 Mg 

4 Mg Y=0.24, =0.864 Mg 

I I 





Figure 3. He- and H-burning luminosities, and the temperature 
in the He-burning shell as a function of time after the first thermal 
pulse for 3 Mq models at F = 0.24 and 0.35, and the 4 Mq, 
Y = 0.24 model. 


(with rju = 0.02), which steeply increases with luminosity. For 
this reason, while the 3 and 4 Mq models of VD09 experience 
a greater number of thermal pulses than our models, the 
opposite is the case with 5 Mq and higher masses, for which 
our models experience a greater number of thermal pulses. 

For the 3 Mq models, variation to the initial He abun¬ 
dance has a small (less than ten per cent) effect on the total 
number of thermal pulses for Y — 0.24, 0.30, and 0.35. How¬ 
ever, with Y = 0.40, the number of thermal pulses almost 
doubles to 48 in comparison with 26 for Y = 0.35. The large 
increase in the number of thermal pulses between Y — 0.35 
and Y = 0.40 is due to the AGB lifetime shortening by 
about 15 per cent, while the average interpulse time drops 
by roughly 50 per cent. 

We use the standard dredge-up efficiency parameter de¬ 
fined by A = AMdredge/AMc, where AMdreUge is the mass 
of He-intershell material mixed to the surface by TDU and 



Figure 4. The maximum He-burning luminosity versus the H- 
exhaused core mass for each thermal pulse of the models at 3, 4, 5, 
and 6 Mq (from left to right) with He mass fractions of T = 0.24 
(black squares), 0.30 (green diamonds), 0.35 (pink circles), and 
0.40 (orange triangles). 


AMc is the growth of the H-exhausted core during the pre¬ 
ceding interpulse phase. 

The TDU is generally less efficient after weaker ther¬ 
mal pulses, and models with a higher core mass attain a 
lower peak He-flash luminosity ( Sackmann|| 1980a Boothroyd 


fc Sackmann|[l988[ [Straniero et al.||2003 1. This is true for 

changes to either the initial mass or the initial He mass 
fraction (Figure]^. However, when models of the same core 
mass are compared, the peak He-burning luminosity in the 
He-enhanced models is less than the primordial-He models 
by up to about 80 per cent. As a consequence, the TDU is 
less efficient in He-enhanced models. 

Figure|^shows the total mass dredged-up by TDU (here¬ 
after TDU mass; equal to the sum of AMdredge) during the 
AGB phase as a function of initial mass for several initial He 
abundances. Increasing the initial He mass fraction typically 
reduces the TDU mass for initial masses from 3 to 6 Mq, 
similar to what KMN14 found with their He-enhanced 1.7 
Mq and 2.36 Mq models. For Y = 0.24, 0.30, and 0.35, the 
maximum dredge-up mass is obtained with an initial mass 
of around 4 Mq. With Y = 0.40 however, the distribution 
changes shape with a peak shifting downwards to approxi¬ 
mately 3 Mq. In contrast with our results, the SCP14 models 
predict that the TDU mass decreases monotonically with 
initial mass, and is overall much lower than for our Y = 0.24 
models (by a factor of up to seven at 4 Mq). The large 
variations between results of SCP14 and the present work 
reflect the significant modelling uncertainties that affect the 
efficiency of the third dredge-up (for a comparison of TDU 
efficiency between different evolution codes, see |Lugaro et al.| 
|2003[ ). Our specific predictions for the TDU efficiency are 
influenced by the mass and time resolution of the extremely 
He-enhanced models, and our treatment of the border be¬ 
tween radiative and convective zones (as demonstrated by 
Frost fc Lattanzio|1996 1. 
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Figure 5. The total third dredge-up mass during the AGB as a 
function of initial mass with He mass fractions oiY = 0.24 (black 
squares), 0.30 (green diamonds), 0.35 (pink circles), and 0.40 
(orange triangles). Results of the 1.7 Mq and 2.36 Mq models 
are from KMN14. 


The TDU mass of the 3 Mq, Y = 0.40 model is unusual 
because it exceeds the TDU mass of the 3 Mq model with 
Y = 0.35 despite the increased He abundance, which reduces 
the dredge-up efficiency. Indeed, both the average and maxi¬ 
mum values of A are lower in the 3 Mq model with Y = 0.40 
than Y = 0.35 (Table [^. The unusually high TDU mass 
with Y = 0.40 is the result of a near-doubling in the number 
of thermal pulses, which leads to more thermal pulses with 
TDU. 

3.4 Hot bottom burning 

At Z = 0.0006 and primordial He abundance, AGB stars 
with initial masses > 4 Mq experience sufficiently high 
temperatures at the base of the convective envelope (Tbce) 
for proton-capture nucleosynthesis to take place there. This 
is called hot bottom burning (HBB). 

With temperatures in the envelope above about 50 MK, 
the resulting activation of the CN cycle begins to convert 
a significant fraction of C nuclei into N. In the absence of 
other effects, the CN-cycle would cause the C/0 ratio in the 
envelope to decrease. However, additional primary nuclei 
from the He-burning shell are periodically transported to 
the envelope via TDU, leading to a surface C/O ratio that 
depends on the interplay between TDU and HBB. In the 
more massive AGB stars at solar metallicity, HBB is obser- 
vationally confirmed to prevent the surface abundances from 
becoming C-rich (C/O > l;|Lattanzio|1992||Boothroyd et al. 


HBB in nature is confirmed by observations of C-deficient, 
N-rich AGB stars in the Magellanic clouds (|McSaveney et al.| 
2007t . 

Figure]^ shows the maximum temperature at the base 
of the convective envelope in our models as a function of 
the initial mass. The higher temperatures in models with 


19931. Further, the combined operation of both TDU and 


Figure 6. The maximum temperature at the base of the convec¬ 
tive envelope during the interpulse phase as a function of initial 
mass with He mass fractions of U = 0.24 (black squares), 0.30 
(green diamonds), 0.35 (pink circles), and 0.40 (orange triangles). 
Results of the 1.7 Mq and 2.36 Mq models are from KMN14. 
The horizontal dashed line at 50 MK indicates the approximate 
temperature above which hot bottom burning significantly alters 
surface abundances. 


increased initial He abundance enables HBB to take place 
at lower initial masses than at primordial He abundance. At 
masses above about 4 Mq, HBB will proceed at increased 
rates with He-enhancement, while models of less than about 
2 Mq will not experience HBB temperatnres even with He 
abundances as high as U = 0.40. For the 3 Mq models, 
Tbce/ increases from 43 MK with Y = 0.24 to 86 MK with 
Y = 0.40. This corresponds to a change from virtually no 
HBB to significant CN cycling in the envelope. 

The maximum envelope temperatures of the primordial- 
He models are in reasonable agreement (less than ten per 
cent difference) with the models of VD09 for masses of 4 
Mq and above. The difference is more significant at 3 Mq, 
where our model has a maximum temperature of 57 MK 
during TDU and 43 MK during the interpulse phase, while 
the corresponding model of VD09 has a maximum temper¬ 
ature of 75 MK. The higher envelope temperatures in the 
models of VD09 are a result of the FST convection model, 
which predicts more efficient convective transport than the 
MLT theory. In contrast to the models of VD09 and the 
models presented here, the models of SCP14 predict much 
lower temperatures at the base of the convective envelope 
(22 MK at 4 Mq versus 84 MK in our model and 91 MK 
for VD09). SCP14 compare their models to the models of 


D’Orazi et al. (20131 (calculated with the same evolutionary 


code used here), which experience more efficient HBB. They 
attribute the difference to some combination of the equation 
of state, interpolation of the radiative opacity tables, and 
their particular mixing scheme. 

Figure]^ shows the maximum temperature at the base 
of the convective envelope and the surface C/O ratio for each 
thermal pulse of the 3 Mq models. Over the first ten or so 
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dances of other light elements through the NeNa cycle and 
the MgAl chain. Activation of the NeNa cycle can produce 



Figure 7. The maximum temperature at the base of the convective 
envelope during the preceding interpulse phase (top panel) and 
the surface C/O ratio (bottom panel) as a function of thermal 
pulse number for 3 Mq models with He mass fractions oiY = 0.24 
(black squares), 0.30 (green diamonds), 0.35 (pink circles), and 0.40 
(orange triangles). The dashed line in the lower panel indicates a 
C/O ratio of 1. 


TPs, the shallower initial rise of C/O in the He-enhanced 
models can be attributed to the less efficient TDU, as none 
of the models have reached HBB temperatures at this stage. 
In models with He-enhancement, HBB becomes active and 
is particularly noticeable with Y = 0.35 and Y = 0.40, for 
which significant CNO cycling in the envelope causes the 
C/O ratio to begin decreasing at around TP 16. As mass loss 
erodes the envelope to below a critical mass value (around 
1.5 Mq in this case) over the last few thermal pulses, HBB 
ceases while dredge-up continues to take place, causing an 


upturn in the C/O ratio in these two models (see Frost et al. 
19981. In the Y = 0.40 model, the cessation of HBB causes 


the surface to transition from O-rich to C-rich (for a second 
time) at the third-last thermal pulse. 

Aside from CNO abundances, HBB can alter the abun¬ 


or destroy ^Na at the expense of Ne isotopes (Arnould et al 


1999 Mowlavi| 1999b Karakas &; Lattanzio|20^ l. At tem¬ 

peratures of 50-200 MK, proton captures on to ^'^Na produce 
more ^°Ne than ^^Mg, as the rate of the ^^Na(p, o l^^Ne reac- 
tion is several times faster than ^®Na(p, 7 )^‘*Mg (Hale et al. 


20041. Activation of the MgAl chain can result in a net pro 
duction of ^®A1 and ^^Al (Arnould et al.|l999 Denissenkov & 


Herwig|2003 Ventura fc D’Antona|2008| ), and further proton 

capture in the most massive AGB stars can produce ^®Si 
(Ventura et al.|2011l. 


3.5 Core mass-luminosity relation 

A well-known correlation exists between the surface lumi¬ 
nosity of an AGB star and its H-exhausted core mass. A 


linear form was first proposed by Paczyriski (19701, although 


stars undergoing HBB and TDU were found to diverge from 
this relation (Bloecker &: Schoenberner|1991 Lattanzio|1992 


Boothroyd fc Sackmann||1992 Marigo et al.||199'9 ). Figure 

^ shows the surface luminosity versus core-mass at each 
thermal pulse for the models of 3, 4, and 5 Mq. 

We also compare the core-mass-luminosity behaviour 
of our models with values from the fitting formula specified 


by equations 29 to 34 of Izzard et al. (20041. This formula 


estimates the surface luminosity as a function of the core 
mass, the envelope mass, and the total growth in the core 
mass during the TP-AGB (neglecting decreases from TDU). 
The formula is based on the stellar evolutionary models of 
[Karakas et al.| ( |2002| ) with He mass fractions that range 
from 0.24 and 0.30, depending on the metallicity. The grid 
of evolutionary models used to generate the fit extends down 
to Z — 0.0001, but only with masses of 1.25 to 2.25 Mq. 
At ^ = 0.004 and higher metallicities, the grid includes 
models with initial masses up to 6 Mq . Although the updated 
formula of jlzzard et al.j ( |2006P improves the fit to models of 
low metallicity, its parameterisation in terms of the initial 
mass causes it to be less accurate for He-rich models, which 
have different stellar structures compared to primordial-He 
models at a given initial mass. 

For the 3 Mq models, the formula is a good fit for 
He mass fractions between Y = 0.24 and extrapolation up 
to y = 0.35. However, our 3 Mq model with Y = 0.40 
is consistently more luminous (by up to 25 per cent) than 
the formula value. Between Y = 0.35 and Y = 0.40, the 
maximum temperature at the base of the convective envelope 
during H-burning increases from 63 to 86 MK, and the more 
extreme HBB further increases the surface luminosity. The 
excess envelope luminosity with Y = 0.40 in comparison 
to the fitting formula (which accounts for the larger core 
mass) is probably partly due to the more He-rich envelope 
composition (which affects the H-burning rate) and also 
because of the smaller envelope mass compared to higher- 
mass models with primordial He and the same core mass. 

The 4 Mq models at all He mass abundances have 
higher peak luminosities than the formula values, although 
this is a small effect with Y — 0.24. The He-enhanced 4 Mq 
models diverge from the relation in a manner characteristic 
of HBB, with a rise in the middle followed by a decline 
as the decreasing envelope mass causes a reduction in the 
temperature at the base of the envelope. A reasonably good 
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Figure 8. Luminosity versus H-exhausted core mass during the 
AGB for 3, 4, and 5 Mq models with He mass fractions of y = 0.24 
(black squares), 0.30 (green diamonds), 0.35 (pink circles), and 
0.40 (orange triangles). The plot also shows values from the fitting 
formula of|Izzard et al.|(|2004^ (1+04). 


fit is found for the 5 Mq models, and this is probably because 
HBB is already active for this initial mass with Y = 0.24. 

In summary, the initial He abundance alters the core- 
mass-luminosity relation, and a more accurate fitting formula 
for application to He-rich populations would require a fit to 
stellar evolutionary models at the appropriate He abundance. 


3.6 Carbon burning and final fates 


The C-burning behaviour and the stellar remnants of the 
models depend on the initial stellar mass and the initial 
He abundance. There is also a dependence on the initial 
metallicity ( Cassisi fc Castellani||19Si3 Umeda et al.|1999 |. 

All of the models presented in this work would eventually 
form white dwarfs (WD) rather than exploding as electron- 
capture supernovae because they end with core masses well 

In 


below the limiting mass of 1.37 Mq (Miyaji et al.| 
Tablej^we classify the stellar remnants for each of our models 
based on the core composition at the end of our calculations, 
which would be similar to the composition of the resulting 
WD remnant. 

AGB stars with initial masses below the minimum mass 
required for off-centre core C ignition (Mup, as defined by 
Becker & Iben 19791 end the AGB with cores composed 


largely of C and O in roughly equal proportions that are set 
by earlier core and shell He burning. Thus, we classify the 
remnants from our models without C ignition as GO WDs. 

In stars with degenerate GO cores following central He 
exhaustion, the position of maximum temperature moves 
outwards from the centre to where the energy liberated by 
gravitational contraction exceeds the rate of cooling by neu¬ 
trino emission ( Becker fc Iben||1979 |. In stars more massive 
than Mup, the temperature maximum exceeds about 600 to 
700 MK ( |Siess|[^06[ ) and G is ignited off-centre in the core. 
In this case, the star is referred to as a ‘super-AGB’ star 
( Garcia-Berro fc Iben||19M Siess||2007| 20101. 

C burning proceeds via the a)^°Ne and 

i2c(12g p)23Na reactions. The latter reaction is followed 
by ^®Na(p, a)^°Ne, and thus also contributes to the abun¬ 
dance of ^°Ne in the core, as does the ^®0(a, 7 )^°Ne reaction. 
The region affected by G burning will be composed largely 
of O and Ne, with a ratio that is affected by the rate of 
7 )^°Ne. The G burning can either be aborted before 


reaching the centre and form a hybrid CO(Ne) WD (Doherty 


et al. poTo] [Karakas et al.p012| [Denissenkov et al.|2013[ ), or 

proceed all the way to the centre and form an ONe WD. 


Bono et al.l (20001 show for super metal-rich compo¬ 


sitions {Z = 0.04) that the Mup boundary decreases with 
increasing He abundance from 9.5 ± 0.5 Mq for T = 0.29 to 
7.7 ± 0.2 Mq for Y — 0.37. With primordial He abundance, 
C burning at Z = 0.0006 takes place with a minimum ini¬ 


tial mass of 6.5 Mq (Doherty et al. 


2015), and we do not 


find G burning in our Y = 0.24 modeis with masses up to 
6 Mq. These models lead to CO WDs. The 6 Mq model 
with Y = 0.30 experiences off-centre C burning that leaves 
C in the inner 0.5 Mq unburned, and this region remains 
composed of C and O. We classify the remnant of this model 
as a hybrid CO(Ne) white dwarf ( Doherty et al.|20T5 ). The 
models at 5 Mq and below with Y = 0.30 do not show any 
significant C burning. 

For Y = 0.35, the 6 Mq model experiences off-centre C 
ignition that proceeds to full central C burning and forms 
an ONe core. Figure shows the C-burning luminosity and 
the behaviour of the C-burning convective zones in this 
model as a function of time. For Y = 0.35, the models at 5 
Mq and below do not show any significant C burning. For 
Y = 0.40, the 5 Mq model experiences off-centre C burning 
that produces a CO(Ne) white dwarf and the 6 Mq model 
undergoes full central C burning that produces an ONe white 
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Figure 9. C-burning characteristics of the 6 Mq, Y = 0.35 model 
versus time before the first thermal pulse. Top: A Kippenhahn 
diagram showing C-burning and envelope convective zones and 
the location of maximum temperature (red). Middle: The total 
luminosity (Ltot) and the luminosities due to He- and C-burning 
(Ljje and Lq). Bottom: The maximum temperature in the model. 


dwarf. The models with Y = 0.40 at 4 Mq and below do 
not show any significant C burning. 

In summary, we find that Mup ed, Z = 0.0006 exceeds 
6 Mq for Y = 0.24, decreases to between 5 and 6 Mq for 
Y = 0.30 to 0.35, and decreases further to between 4 and 
5 Mq with Y = 0.40. Future investigations with a finer 
mass grid would be required to determine more precisely the 
dependence of Mup on the He abundance. 



Miui [Mq] 

Figure 10. Average He mass fractions in the ejecta of the stellar 
models versus initial mass with He mass fractions ol Y = 0.24 
(black squares), 0.30 (green diamonds), 0.35 (pink circles), and 
0.40 (orange triangles). 1.7 Mq and 2.36 Mq models are from 
KMN14. 

4 NUCLEOSYNTHESIS AND STELLAR 
YIELDS 

The stellar yield represents the contribution of a particular 
chemical species i (e.g., an element or nuclide) to the inter¬ 
stellar medium by a star over its lifetime, and is calculated 
using the formula 

[ Xi{t)M{t)dt, (1) 

Jo 

where Xi{t) is the mass fraction of species i at time t and 
M{t) is the stellar mass-loss rate at time t. For models with 
a non-zero envelope mass at the end of our calculations, we 
assume that all mass exterior to the core is ejected with the 
composition of the surface at the last computed model. 

The stellar yields of all models are provided as online 
data tables, with an example of their form and content in 
Table |A1| In the online data tables, we also include the 
net yields, which have the initial abundances subtracted 
according to 

^netyield^ / [Yi (t) - A, (0)] M(t) df, (2) 

Jo 

where Yi(0) is the mass fraction of species i in the initial 
composition. The net yields indicate whether a chemical 
species is produced or depleted in the context of chemical 
evolution. 


4.1 Yields of light elements 

In AGB stars, the primary He produced by H-burning reac¬ 
tions is mixed to the surface during dredge-up events, with 
the largest increase to the surface abundance of He occurring 
during SDU and a smaller increase due to TDU. As we show 
in Figure the average He mass fractions in the ejecta of 
our models are typically 0.05 to 0.10 higher than the initial 
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He mass fractions, with the higher-mass models releasing 
material that is more He-enhanced due to a deeper SDU 
( Karakas et al.|2006 l. 

For primordial-He models, the He mass fractions of our 
yields are generally higher than the predictions of SCP14 or 
VD09. At the largest difference, our predictions for (F^jecta) 
are higher than VD09 and SCP14 by about 0.05 and 0.03, 
respectively. The FRUITY database ( Cristallo et al.|2011 | 
provides surface He mass fractions of the SCP14 models 
after each dredge-up event, which enables us to compare the 
effects of SDU and TDU. At 3 Mq, the surface Y following 
SDU is 0.30 (cf., 0.26, SCP14) and increases with each TDU 
episode up to 0.33 at the end of our calculations, while the 
model of SCP14 finishes with a surface Y of 0.27. The higher 
He mass fraction in our yield at 3 Mq compared with SCP14 
is caused by a deeper SDU and more extensive TDU in our 
model. For initial masses of 4 to 6 Mq, our post-SDU He 
abundances are very similar to the models of SCP14, and 
the differences in (Tejecta) largely reflect differences in the 
efficiency of HBB. 

Figure ED illustrates the time evolution of the surface 
abundances of Ne, Na, Mg, Al, and Si isotopes in the 4 
Mq models with Y = 0.24 (left panel) and Y = 0.40 (right 
panel). We focus on the 4 Mq models as they demonstrate 
the significant changes that can occur due to HBB with 
variations to the initial He abundance. As shown in Table 
the maximum temperature at the base of the envelope during 
the interpulse phase increases from 84 MK at Y = 0.24 to 101 
MK at y = 0.40, with implications for the nucleosynthesis 
via proton capture reactions. However, while the maximum 
temperature increases, other important factors include the 
length of time during which the envelope is subject to high 
temperatures and the total mass of TDU (tnBB and M^u 
in Table [^, both of which are lower in the Y = 0.40 model. 


The surface abundance of ^^Ne displayed in Figure 
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for both y = 0.24 and Y = 0.40 follows a stair-step increase 
as ^^Ne is periodically dredged up into the envelope following 
its production in convective pulses. The surface abundance of 
^^Na, which increases monotonically at Y = 0.24, is initially 
destroyed by proton captures in the envelope at Y = 0.40. 
Nucleosynthesis in the convective pulses also contributes 
to the intershell abundances of ^®Mg and ^®Mg, which are 
produced via the ^^Ne(Q;, n)^®Mg and ^^Ne(a, 7 )^®Mg reac¬ 
tions. At Y = 0.24, increases to ^®Mg and ^®Mg occur during 
dredge-up events, with minimal increases due to HBB. At 
Y = 0.40, however, proton captures cause the ^®Mg abun¬ 
dance to increase rapidly at the expense of ^^Mg. While some 
of the ^“^Mg is destroyed by proton captures at Y = 0.24 
(0.1 dex in mass fraction), the destruction of ^"^Mg is far 
more extensive at Y = 0.40, with a decrease by almost 3.0 
dex before the reduced envelope temperatures and efficient 
dredge-up cause the abundance to begin increasing. The 
production of ^^Al is lower at Y = 0.40 despite the higher 
TbcT, which indicates the greater impact of a shorter AGB 
phase and reduced dredge-up in the He-rich model for this 
nuclide. In both cases, the amount of ^®Si produced from 
proton capture on to ^^Al is very small (< 0.1 dex). 

Figure 12 shows the stellar yields (in Mq) of the light 


around 3 Mq. This is the highest mass model that does not 
experience HBB, while HBB is highly effective in reducing 
the C yield and increasing the yield of N with masses of 
about 4 Mq and above. With helium mass fractions from 
0.30 to 0.40, the model that produces the highest C yield 
becomes the 2.36 Mq model of KMN14. For the models of 3 
Mq, the increase in the N yield with He abundance shows 
that the increased efficiency of HBB with He-enhancement 
overcomes the opposing effect of reduced TDU. However, 
this changes with M > 4 Mq , where the yields of N decrease 
with increasing He, demonstrating that the reduced TDU 
and shorter AGB phases have a greater effect on the N yields 
than the increased temperatures at the base of the envelope. 

For most of the models, the sum of the G-l-N-l-O yields 
is substantially reduced by helium enhancement, with the 
exception of the 6 Mq model with Y = 0.40. The 6 Mq 
model with Y = 0.40 experiences a corrosive SDU that 
dredges up C and O from the core, which affects the stellar 
yields of several light elements. Compared with the 6 Mq 
model at Y = 0.35 (which does not experience a corrosive 
SDU), the yield of N (produced from dredged-up C) increases 
by 70 per cent, while O increases by 1600 per cent and the 
Mg yield increases by 80 per cent. 

The predicted yields of C for the models of SCP14 are 
up to three times lower than our Y = 0.24 models at masses 
of 3 Mq and below, and roughly two times higher for masses 
of 4 Mq and above. The smaller C yields in the lower mass 
range reflect the significantly reduced TDU masses of the 
SCP14 models, which are up to 85 per cent lower than the 
models presented here. The higher C yields of the SCP14 
models in the upper mass range occur in spite of the reduced 
TDU and are caused by lower envelope temperatures and 
consequently less HBB than the models presented here. 

The only stable isotope of F is ^®F, which is produced 
via ^'*N(q:, 7)^®F((5+z/)^®0(p, a)^®N(Q;, 7 )^®F during thermal 
pulses in the He-intershell ( |Jorissen et al.|1992 1. F can be de¬ 
stroyed both by proton captures in the envelopes of stars with 
HBB via ^®F(p, a)^®0 and a captures in the He-intershells 
of stars more massive than about 5 Mq via ^®F(a,p)®®Ne 
( Lugaro et al.|20M Gristallo et al.|2014 l. Thus, the F yields 
depend on both the TDU mass and the temperatures at the 
base of the convective envelope, causing a narrow peak of 
production at around 2 to 3 Mq that makes F abundances a 
powerful indicator of stellar mass (e.g., D’Orazi et al.|2013 l. 

The yield of Ne is sensitive to the inclusion and extent of 
a PMZ ( Shingles &: Karakas|2013 1. This is because additional 
primary ^'‘N produced in the PMZ is converted into Ne during 
convective pulses via ^‘*N(a, 7)^®F(/I"'"j/)^®0(a, 7 )®®Ne. The 
very short half-life of ^®F (110 min) means that the additional 
F nuclei produced through this chain will not reach the 
surface (unless neutrons are captured to make ^®F). The 
increase in the yield of F with the inclusion of a PMZ is caused 
by the additional production of ^®N, which is converted into 
^®F during ther mal pulses. 

shows the isotopic fractions of ®‘^Mg, 
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"Mg, 


Figure 

and ®®Mg in the stellar yields as a function of the initial 
stellar mass and He abundance. The shaded regions indicate 
the approximate range of Mg isotopic fractions for red giants 


et al. (20131. 


elements C^N, O, F, Ne, Na, Mg, and Al as a function of in ui Centauri with [Fe/H] > —1.4, as measured by Da Costa 
initial mass. 

At primordial He abundance, the yield of C increases 
with the initial mass (and TDU mass) up to a peak at 


‘Mg 


is efficiently produced by the a-process in core-collapse 
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Figure 11. Surface abundance of selected nuclides that participate in the NeNa cycle and the MgAl chain versus time in 4 Mq models 
at y = 0.24 (left panel) and Y = 0.40 (right panel). Time on the horizontal axis is relative to the time of the first thermal pulse. 


supernovae, while the presence of and ^®Mg at low 

metallicities is generally a good indicator of enrichment by 
massive AGB stars or Wolf-Rayet stars. In the u Centauri 
stars, the range of ^'*Mg/Mg fractions of about 50 to 80 
per cent is lower than would be expected from a pure core¬ 
collapse enrichment scenario (e.g., > 90 per cent for the 
low-metallicity yields of Kobayashi et al.|2011 1 and includes 
values well below the solar value of 79 per cent ([Asplund et al.| 


20091. However, a precise accounting for the contribution by 
AGB stars to these low ^'*Mg/Mg values in uj Centauri would 
require a detailed chemical evolution model (e.g., [Romano| 


et al. 20071, and it is not immediately clear whether the 


yields of our He-rich models provide a better match to these 
observations. 

This use of our yields or those of VD09 could potentially 
constrain the mass range of AGB polluters in to Centauri. 
This is because the Mg isotope ratios in these yields vary 
strongly as a function of the initial stellar mass. In contrast, 
the limited HBB in the models of SCP14 leads to low isotopic 
fractions of ^®Mg and ^®Mg, which remain similar to the 
values in the initial composition throughout the mass range 
from 1 to 6 Mq. These models predict that Mg isotope ratios 
provide only a very weak constraint on the mass range of 
AGB polluters. 

In our models for Y = 0.40, the ^®Mg/Mg ratio sharply 
decreases with initial masses greater than about 3 Mq. The 
envelope temperatures in these models enter the temperature 
range (> 110 MK) where the rate of the ^®Mg(p, 7)^^A1 
reaction begins to exceed that of ^®Mg(p, 7)^®A1. A similar 
trend is not seen at lower He abundances. 


4.2 Yields of neutron-capture elements 

Figurej^shows the heavy element abundance distributions in 
the average of the ejecta for the 3 Mq and 6 Mq models. The 
initial masses of 3 and 6 Mq models are chosen to represent 


the nucleosynthetic behaviour of low- and intermediate-mass 
AGB models. 

In the 3 Mq models, the overall shape of the abundance 
distribution is relatively unchanged by variation to the initial 
He abundance. The primordial model is in reasonable agree¬ 
ment with the shape of the distribution predicted by SCP14, 
but we predict higher absolute abundances at the first and 
second s-process peaks near Y and Ba by around 0.5 dex 
(as represented by [X/Fe] ratios). The greater abundances 
of these elements in our model compared with SGP14 is 
caused by a significantly higher TDU mass (by a factor of 
four) and possibly also our different approaches for modelling 
^®G-pocket formation. 

In contrast with the 1.7 Mq and 2.36 Mq yields of 
KMN14 that steadily decrease with increasing He mass frac¬ 
tion, the Ba, La, Ge, and Pb yields of our 3 Mq models 
exhibit a complicated dependence on the initial He content. 
In the 3 Mq model with Y = 0.40, the heavy-element abun¬ 
dances near the second s-process peak at Ba are higher than 
the model with Y = 0.24, despite the lower TDU mass. This 
is probably a consequence of the greater number of thermal 
pulses in the Y = 0.40 model (48 versus 26 with Y = 0.24) 
and a more significant effect of the ^^Ne source. 

The maximum extent of the pulse-driven convective 
zone (M™(f,l). in Table gives an indication of the size of 
the He-rich intershell, which decreases with increasing He 
content. For the 3 Mq models, the PMZ mass of 10“® Mq 
increases as a fraction of the He-intershell mass from 10 per 
cent with Y = 0.24 to 33 per cent for Y = 0.40. However, 
as discussed in Section the high temperatures at the 
base of the envelope probably prevent a significant s-process 
production from pockets in 3 Mq stars with Y = 0.40. 
We expect that the model for Y — 0.40 with no PMZ is more 
realistic, although there could be a very small contribution 
from pockets in stars of this mass and He abundance. 

For the 6 Mq models, the reduction in TDU mass with 
He-enhancement causes heavy-element abundances in the 
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Figure 12. Yields of selected light elements as a function of initial mass with He mass fractions of Y = 0.24 (black squares), 0.30 (green 
diamonds), 0.35 (pink circles), and 0.40 (orange triangles). Open points represent 3 Mq models without a PMZ. Models with masses of 3 
Mq and below otherwise include a PMZ of 1 X 10~^ 2.36 Mq yields are from KMN14. 


yields to decrease monotonically with increases to the He 
abundance. In models with M > 4 Mq, the main neutron 
source is the ^^Ne(a, n)^®Mg reaction operating in convective 
pulses, which produces much lower neutron-to-Fe-seed ratios 
and higher neutron densities compared to the neutron 
source. As a result, the 6 Mq abundance distributions peak 
at Rb, with positive [Rb/Sr] and [Rb/Zr] ratios that are 
characteristic of s-process production via the neutron 
source. 


This is because the abundance ratios [Rb/Zr] and 
[Rb/Sr] depend on the s-process branchings at ^^Kr (half- 


life = 10.8 yr) and ®®Rb (half-life = 18.6 days) (van Raai 
|et al.|[2012| [Karakas et al.|[2012p . With neutron densities 
above 10“-10“ cm“'^, neutron captures on to ®®Kr and ®®Rb 


cause the s-process path to produce ®^Rb. This nuclide has 
a magic number of neutrons, which gives it a very small 
neutron-capture cross section relative to ®®Rb and neigh¬ 
bouring nuclides and hence it accumulates ( Heil et al.|2008 1. 
At high neutron densities, this causes a larger production 
of Rb than Sr and Zr, which are not affected by branching 
points (D’Orazi et al.|2013|. 


Pb and Bi represent the end point of the s-process chain, 
because heavier nuclei are unstable against a- and /3-decays. 
Pb is not produced in significant quantities by the models 
that do not have a PMZ, as the low neutron-to-Fe seed ratios 
in these models are not sufficient to populate the s-process 
chain to its heaviest nuclide. 


Figure [T^ shows the stellar yields of selected s-process 
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Figure 14. Average abundances in the ejecta versus atomic number for the 3 Mq models (top panel) and 6 Mq models (bottom panel) 
for He mass fractions of Y = 0.24 (black squares), 0.30 (green diamonds), 0.35 (pink circles), and 0.40 (orange triangles). Open points 
indicate 3 Mq models without a PMZ. Models with initial masses of 3 Mq and below otherwise include a PMZ of 1 X 10“® ^^Q- 


elements as a function of initial mass for all of our models. 
Overall, the yields of s-process elements exhibit a very strong 
reduction with He-enhancement, mostly because of the lower 
TDU mass. 

The yields of Ba, La, Ce, and Pb increase with initial 
mass up to a peak near 3 Mq. The initial rise is due to the 
correlation between initial mass and the amount of TDU 
during the AGB. The trend then reverses beyond 3 Mq due 
to our assumption of no PMZs (and pockets) in the 
higher-mass models. 

For our models between 4 and 6 Mq, in which the main 
neutron source is the ^^Ne(Q, n)^®Mg reaction operating in 
convective pulses, the elements Ba, La, Ce, and Pb are not 
produced in significant quantities. With He-enhancement, the 
yields of these elements are reduced further to effectively zero. 
For the lighter s-process elements that are mainly produced 
by models in this mass range (Kr, Rb, Sr, Y, and Zr) the 
yields are decreased by an order of magnitude with Y = 0.40. 

In contrast with our results for Y = 0.24, the SCP14 
predictions for s-process yields are substantially lower. This 
is connected with the smaller number of thermal pulses in 
their models, which result in fewer neutron producing events 
and TDU episodes, and possibly the different treatment of 
^®C-pocket formation. The smaller variation of the yields as 


a function of initial mass for the models of SCP14 reflects 
the flatter profile of TDU mass versus the initial mass. 


The top panel of Figure 16 shows the average [Rb/Zr 


ratio in the ejecta as a function of the initial stellar mass. 
The transition from neutrons produced in pockets to 
neutrons produced by the ^^Ne source in convective pulses 
(which produce higher neutron densities) at around 3 Mq 
is marked by [Rb/Zr] values that are negative or positive, 
respectively. In contrast with the stellar yields of Rb and 
Zr (which vary by more than 1 dex between our models 
and SCP14), our predictions for [Rb/Zr] agree with those 
of SCP14 to within 0.2 dex for initial masses between 3 
and 6 Mq. This is because the [Rb/Zr] ratio is relatively 
independent of the TDU mass, which is substantially different 
in the SCP14 models. 

The relative behaviour of the three s-process peaks is 
quantified by the Is and hs indices and the Pb abundances. 


Following Cristallo et al. (20111, we define [Is/Fe] by 


[Is/Fe] = ([Sr/Fe] -k [Y/Fe] -k [Zr/Fe])/3, 


(3) 


and [hs/Fe] by. 


[hs/Fe] = ([Ba/Fe] -k [La/Fe] -k ]Nd/Fe] -k [Sm/Fe])/4. (4) 
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Figure 15. Yields of selected neutron-capture elements as a function of initial mass with He mass fractions of Y = 0.24 (black squares), 
0.30 (green diamonds), 0.35 (pink circles), and 0.40 (orange triangles). Symbols are the same as Figure [T^ 


Figure 16 shows the average [hs/lsQand [Pb/hs] indices 
in the ejecta of our models as a function of the initial stellar 
mass. Both the [hs/ls] and [Pb/hs] indices are highly sensitive 
to the inclusion of a PMZ (and pockets), which results 
in positive values of these ratios for initial masses less than 
about 3 Mq and negative values for higher masses. 

The [hs/ls] and [Pb/hs] indices of the SCP14 models 
decrease more slowly with initial mass than the models pre¬ 
sented here. This is because their treatment of the PMZ does 
not include an explicit cutoff at a particular initial mass. The 
models of SCP14 include a PMZ that gradually decreases 
in size and neutron production as the initial mass increases. 


^ [hs/ls] = [hs/Fe] — [Is/Fe] 


whereas our models with masses above 3 Mq do not include 
a PMZ. 


5 DISCUSSION & CONCLUSIONS 

We have presented and analysed the results of calculations 
of intermediate-mass stellar models at low metallicity with 
initial He mass fractions of Y = 0.24, 0.30, 0.35, and 0.40. 

The results of this study support the conclusion that the 
initial He abundance is a crucially important parameter for 
stellar evolution and chemical yields. For example, we have 
shown that the dredge up efficiency and the total TDU mass 
in intermediate-mass models are significantly reduced by 
He-enhanced initial compositions. In synthetic AGB models 
that require the value of A as a function of mass and metal- 
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Figure 13. The isotopic fractions (by number) of Mg in the 
stellar yields as a function of initial mass for He mass fractions 
of y = 0.24 (black squares), 0.30 (green diamonds), 0.35 (pink 
circles), and 0.40 (orange triangles). Open points indicate 3 Mq 
models without a PMZ. Models with initial masses ^ 3 Mq 
otherwise include a PMZ of 1 x 10“^ ^©- The points at 1.7 
Mq and 2.36 Mq are from the models of KMN14 (which have 
[a/Fe] = 0.4 for Y = 0.24). The grey shaded regions indicate the 
approximate range of red giants in lj Centauri with [Fe/H] > —1.4 
from observations by |Da Costa et al.| | |201^ . 

licity (e.g., [Bertelli et al.|2008[|BuelH2013[ ), an improvement 
in accuracy would be achieved if the He abundance were 
incorporated as a third parameter and values were obtained 
from full stellar evolutionary models with the appropriate 
initial He abundance. 

One of the objectives of this study was to discover 



Figure 16. The s-process indices [Rb/Zr] (top), [hs/ls] (middle) 
and [Pb/hs] (bottom) of the yields as a function of initial mass 
with He mass fractions of Y = 0.24 (black squares), 0.30 (green 
diamonds), 0.35 (pink circles), and 0.40 (orange triangles). Symbols 
are the same as Figure [T^ 

whether the reduction in the yields of Ba-peak elements 
found by KMN14 with He-enhanced 1.7 Mq and 2.36 Mq 
models is also the case for higher initial masses. For 4 Mq 
and higher masses, we find that it is the case that increasing 
the He mass fraction results in lower yields of Ba, La, and 
Ce. However, our models at 4 Mq and above make a small 
contribution to these elements with Y = 0.24, and hence 
the reduction in these yields with He-enhancement will have 
minimal impact on chemical evolution. Of greater impor¬ 
tance are the s-process yields of models with M < 3 Mq, 
which are subject to the major uncertainties of ^®C-pocket 
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formation. If we assume that the same PMZs form after each 
thermal pulse with TDU for Y = 0.40 and Y = 0.24, then 
predicted Ba-peak yields are higher by almost a factor of 
two, due to the near-doubling of the total number of thermal 
pulses in the He-enhanced model. However, the very high 
temperatures in the envelope during TDU and the thinner 
He-intershell (by a factor of three) for Y = 0.40 suggest 
that a future, more advanced treatment of the PMZ would 
result in a much smaller production of free neutrons for the 
s-process. In this case, the production of Ba-peak elements 
and Pb would be significantly lower at high He abundances 
because of the reduced contribution from pockets. 

Although many conclusions can be drawn from a grid of 
individual stellar yields, an understanding of their combined 


(e.g., Travaglio et al.||1999 

Kobayashi et al.|2011 

Bisterzo 

et al. 2014 Shingles et al. 

20141. From the stellar yields 


of the 1.7 Mq and 2.36 Mq models by KMN14 and the 3 
to 6 Mq models presented here, the chemical evolution of 
s-process elements near [Fe/H] « —1.4 could be predicted as 
a function of He abundance, once other chemical evolution 
assumptions have been made (e.g., Romano et al.|2007 for 
light and a elements). We expect that a comparison between 
the chemical evolution predictions with yields oi Y = 0.40 
models and the s-process abundances in w Centauri would 
be particularly insightful. 

We have compared the stellar yields of our models with 
the primordial-He models of SCP14 and VD09, but the 
present lack of stellar yields for He-rich compositions makes 
it difficult to understand how the predicted impact of He- 
enhancement would be different with alternative modelling 
assumptions and stellar evolution codes. Of critical impor¬ 
tance for the stellar yield predictions is the mass of intershell 
material that is transported to the surface by TDU. The find¬ 
ing of KMN14 and this work that the total mass dredged-up 
by TDU is highly reduced by He-enhancement (up to 96 per 
cent at 6 Mq) needs to be investigated with other stellar 
evolution codes, which already predict much less efficient 


dredge-up with primordial He abundance (see e.g., Mowlavi 
1999a Lugaro et al.|[2003 20121. 


The increase of the core mass with initial He content 


is already a well-established prediction (e.g., Becker & Iben 
|1979| [Lrattanzio||1986[ ), and therefore we expect that any pre¬ 
dictions that are a direct consequence of larger core masses 
(e.g., shorter interpulse periods, shorter AGB lifetimes, and 
higher luminosities) would be qualitatively similar with al¬ 
ternative modelling assumptions. However, the predictions 
of AGB lifetimes and the total number of thermal pulses are 
highly dependent on the chosen prescription for the mass- 
loss rate. The increased mass-loss rate in our models with 
He-rich compositions would occur to a different extent with 
alternative mass loss prescriptions, which have a different de¬ 


pendence on the stellar radius and luminosity (e.g., Bloecker 
1995 Straniero et al.|2006 l. 


The changing fates of stars with the same initial mass 
but different He abundances leads to questions about how 
other mass boundaries are shifted by He-enhancement. For ex¬ 
ample, the minimum initial mass to form an electron-capture 
supernova at ^ = 0.0006 with primordial He abundance is 
about 8.2 to 8.4 Mq (Doherty et al.|2015l. With the larger 
core masses of He-enriched stars, this boundary would shift 
to lower masses, which are more numerous with a standard 


initial mass function. As a consequence, the rates of electron- 
capture supernovae and neutron star formation would be 
higher at a given star-formation rate with He-enhancement. 

Similarly, He-enhancement would increase the number of 
hybrid GO(Ne) WDs, which have been suggested as possible 
progenitors of Type lax supernovae (Denissenkov et al.|2015 


[Kromer et al.|20T^ . If this is the case, then the rate of Type 
la supernovae would be higher for He-rich populations. 


ACKNOWLEDGMENTS 


We thank the anonymous referee for helpful corrections to 
the manuscript. LJS thanks Cherie Fishlock for proofreading 
the manuscript. AIK was supported through an Australian 
Research Gouncil Future Fellowship (FTl 10100475). RJS 
is the recipient of a Sofja Kovalevskaja Award from the 
Alexander von Humboldt Foundation. ML is a Momentum 
Project Leader of the Hungarian Academy of Sciences. This 
research was undertaken with the assistance of resources from 
the National Computational Infrastructure (NCI), which is 
supported by the Australian Government. This research 
has made use of NASA’s Astrophysics Data System. The 


MATPLOTLIB package was used to generate plots (Hunter 
200711. 


REFERENCES 

Abia C., Dominguez L, Gallino R., Busso M., Masera S., 
Straniero O., de Laverny P., Plez B., Isern J., 2002, ApJ, 
579, 817 

Arnould M., Goriely S., Jorissen A., 1999, A&A, 347, 572 
Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, 
ARA&A, 47, 481 

Becker S. A., Iben I. J., 1979, ApJ, 232, 831 
Bedin L. R., Piotto G., Anderson J., Cassisi S., King I. R., 
Momany Y., Carraro G., 2004, ApJ, 605, L125 
Bertelli G., Girardi L., Marigo P., Nasi E., 2008, A&A, 484, 
815 

Biermann L., 1948, Zeitschrift fiir Astrophysik, 25, 135 
Bisterzo S., Gallino R., Straniero O., Cristallo S., Kappeler 
F., 2012, MNRAS, 422, 849 

Bisterzo S., Travaglio C., Gallino R., Wiescher M., Kappeler 
F., 2014, ApJ, 787, 10 
Bloecker T., 1995, A&A, 297, 727 
Bloecker T., Schoenberner D., 1991, A&A, 244, L43 
Bonacic Marinovic A., Lugaro M. A., Reyniers M., van 
Winckel H., 2007, A&A, 472, LI 

Bono G., Caputo F., Cassisi S., Marconi M., Piersanti L., 
Tornambe A., 2000, ApJ, 543, 955 
Boothroyd A. L, Sackmann I.-J., 1988, ApJ, 328, 671 
Boothroyd A. L, Sackmann 1. J., 1992, ApJ, 393, L21 
Boothroyd A. L, Sackmann I.-J., Ahern S. C., 1993, ApJ, 
416, 762 

Buell J. F., 2013, MNRAS, 428, 2577 
Busso M., Gallino R., Wasserburg G. J., 1999, ARA&A, 37, 
239 

Campbell S. W., D’Orazi V., Yong D., Constantino T. N., 
Lattanzio J. C., Stancliffe R. J., Angelou G. C., Wylie-de 
Boer E. C., Grundahl E., 2013, Nature, 498, 198 
Cannon R. C., 1993, MNRAS, 263, 817 







































He-rich AGB evolution and nucleosynthesis 19 


Canute V. M., Mazzitelli I., 1991, ApJ, 370, 295 
Cassisi S., Castellani V., 1993, A&AS, 88, 509 
Charbonnel C., Chantereau W., Decressin T., Meynet G., 
Schaerer D., 2013, A&A, 557, L17 
Chiappini C., Matteucci F., Romano D., 2001, ApJ, 554, 
1044 

Chiefs A., Straniero O., 1989, A&AS, 71, 47 
Christy-Sackmann I. J., Paezyhski B., 1975, Memoires of 
the Societe Royale des Sciences de Leege, 8, 335 
Cristallo S., Imbriani C., Piersanti L., Abia C., Cialanella 
L., Straniero O., 2014, A&A, 570, 46 
Cristallo S., Piersanti L., Straniero O., Gallino R., 
Dominguez L, Abia C., Di Rico G., Quintini M., Bisterzo 
S., 2011, ApJS, 197, 17 

Cristallo S., Straniero O., Gallino R., Piersanti L., 
Dominguez I., Lederer M. T., 2009, ApJ, 696, 797 
Da Gosta G. S., Marino A. F., 2011, PASA, 28, 28 
Da Gosta G. S., Norris J. E., Yong D., 2013, ApJ, 769, 8 
De Smedt K., van Winckel H., Karakas A. I., Siess L., 
Goriely S., Wood P. R., 2012, A&A, 541, 67 
Denissenkov P. A., Herwig F., 2003, ApJ, 590, L99 
Denissenkov P. A., Herwig F., Truran J. W., Paxton B., 
2013, ApJ, 772, 37 

Denissenkov P. A., Tout G. A., 2003, MNRAS, 340, 722 
Denissenkov P. A., Truran J. W., Herwig F., Jones S., 
Paxton B., Nomoto K., Suzuki T., Toki H., 2015, MNRAS, 
447, 2696 

Doherty C. L., Gil-Pons P., Lau H. H. B., Lattanzio J. C., 
Siess L., Gampbell S. W., 2014, MNRAS, 441, 582 
Doherty C. L., Gil-Pons P., Siess L., Lattanzio J. C., Lau 
H. H. B., 2015, MNRAS, 446, 2599 
Doherty G. L., Siess L., Lattanzio J. G., Gil-Pons P., 2010, 
MNRAS, 401, 1453 

D’Orazi V., Lucatello S., Lugaro M. A., Gratton R. G., 
Angelou G., Bragaglia A., Carretta E., Alves-Brito A., 
Ivans 1. L, Masseron T., Mucciarelli A., 2013, ApJ, 763, 
22 

D’Orazi V., Lugaro M. A., Gampbell S. W., Bragaglia A., 
Garretta E., Gratton R. G., Lucatello S., D’Antona F., 
2013, ApJ, 776, 59 

Dupree A. K., Avrett E. H., 2013, ApJL, 773, L28 
Fishlock C. K., Karakas A. L, Lugaro M. A., Yong D., 2014, 
ApJ, 797, 44 

Fishlock C. K., Karakas A. L, Stancliffe R. J., 2014, MNRAS, 
438, 1741 

Frost C. A., Cannon R. C., Lattanzio J. C., Wood P. R., 
Forestini M., 1998, A&A, 332, L17 
Frost C. A., Lattanzio J. C., 1996, ApJ, 473, 383 
Gallino R., Arlandini C., Busso M., Lugaro M. A., Travaglio 

C. , Straniero O., Chiefii A., Limongi M., 1998, ApJ, 497, 
388 

Garcia-Berro E., Iben L, 1994, ApJ, 434, 306 
Garci'a-Hernandez D. A., Zamora O., Yagiie A., Uttenthaler 

S. , Karakas A. L, Lugaro M. A., Ventura P., Lambert 

D. L., 2013, A&A, 555, L3 

Gil-Pons R, Doherty C. L., Lau H., Campbell S. W., Suda 

T. , Guilani S., Gutierrez J., Lattanzio J. C., 2013, A&A, 
557, 106 

Goriely S., Mowlavi N., 2000, A&A, 362, 599 
Goriely S., Siess L., 2004, A&A, 421, L25 
Hale S. E., Ghampagne A. E., Iliadis G., Hansper V. Y., 
Powell D. C., Blackmon J. G., 2004, Phys. Rev. G, 70, 


045802 

Heil M., Kappeler F., Uberseder E., Gallino R., Bisterzo S., 
2008, Phys. Rev. G, 78, 025802 
Herwig E., 2000, A&A, 360, 952 
Herwig E., 2005, ARA&A, 43, 435 
Herwig E., Langer N., 2001, MmSAI, 72, 277 
Herwig E., Langer N., Lugaro M. A., 2003, ApJ, 593, 1056 
Hollowell D., Iben 1. J., 1989, ApJ, 340, 966 
Hunter J. D., 2007, Computing in Science and Engineering, 

9, 90 

Iben 1. J., 1975, ApJ, 196, 525 
Iben 1. J., Renzini A., 1982, ApJ, 259, L79 
Iglesias C. A., Rogers F. J., 1996, ApJ, 464, 943 
Izzard R. G., Dray L. M., Karakas A. L, Lugaro M. A., 
Tout G. A., 2006, A&A, 460, 565 

Izzard R. G., Glebbeek E., Stancliffe R. J., Pols O. R., 2009, 
A&A, 508, 1359 

Izzard R. G., Tout C. A., Karakas A. L, Pols O. R., 2004, 
MNRAS, 350, 407 

Joo S.-J., Lee Y.-W., 2013, ApJ, 762, 36 
Jorissen A., Smith V. V., Lambert D. L., 1992, A&A, 261, 
164 

Karakas A. L, 2014, MNRAS, 445, 347 
Karakas A. L, Fenner Y., Sills A., Campbell S. W., Lattanzio 
J. C., 2006, ApJ, 652, 1240 

Karakas A. L, Garefa-Hernandez D. A., Lugaro M. A., 2012, 
ApJ, 751, 8 

Karakas A. L, Lattanzio J. C., 2003, PASA, 20, 279 
Karakas A. L, Lattanzio J. C., 2014, PASA, 31, e030 
Karakas A. L, Lattanzio J. C., Pols O. R., 2002, PASA, 19, 
515 

Karakas A. L, Marino A. F., Nataf D. M., 2014, ApJ, 784, 
32 

King 1. R., Bedin L. R., Cassisi S., Milone A. P., Bellini A., 
Piotto G., Anderson J., Pietrinferni A., Cordier D., 2012, 
AJ, 144, 5 

Kobayashi C., Karakas A. L, Umeda H., 2011, MNRAS, 
414, 3231 

Kromer M., Ohlmann S. T., Pakmor R., Ruiter A. J., Hille- 
brandt W., Marquardt K. S., Ropke F. K., Seitenzahl 1. R., 
Sim S. A., Taubenberger S., 2015, MNRAS, 450, 3045 
Langer N., Heger A., Wellstein S., Herwig F., 1999, A&A, 
346, L37 

Lattanzio J. G., 1986, ApJ, 311, 708 

Lattanzio J. G., 1992, Astronomical Society of Australia, 

10, 120 

Lau H. H. B., Gil-Pons P., Doherty C. L., Lattanzio J. C., 
2012, A&A, 542, A1 

Lebzelter T., Wood P. R., 2007, A&A, 475, 643 
Lodders K., 2003, ApJ, 591, 1220 

Lugaro M. A., Herwig F., Lattanzio J. G., Gallino R., 
Straniero O., 2003, ApJ, 586, 1305 
Lugaro M. A., Karakas A. L, Stancliffe R. J., Rijs C., 2012, 
ApJ, 747, 2 

Lugaro M. A., Ugalde C., Karakas A. L, Gorres J., Wiescher 
M., Lattanzio J. C., Gannon R. G., 2004, ApJ, 615, 934 
MeSaveney J. A., Wood P. R., Scholz M., Lattanzio J. C., 
Hinkle K. H., 2007, MNRAS, 378, 1089 
Magic Z., Asplund M., 2014, eprint arXiv:1405.7628 
Marigo P., Aringer B., 2009, A&A, 508, 1539 
Marigo P., Girardi L., Weiss A., Groenewegen M. A. T., 
1999, A&A, 351, 161 



20 Shingles et al. 


Marino A. F., Milone A. P., Przybilla N., Bergemann M., 
Lind K., Asplund M., Cassisi S., Catalan M., Casagrande 

L. , Valcarce A. A. R., Bedin L. R., Cortes C., D’Antona F., 
Jerjen H., Piotto G., Schlesinger K., Zoccali M., Angeloni 

R. , 2014, MNRAS, 437, 1609 

Miglio A., Brogaard K., Stello D., Chaplin W. J., D’Antona 
F., Montalban J., Basu S., Bressan A., Grundahl F., Pin- 
sonneault M., Serenelli A. M., Elsworth Y., Hekker S., 
Kallinger T., Mosser B., Ventura P., Bonanno A., et al. 
2012, MNRAS, 419, 2077 

Miszalski B., Boffin H. M. J., Jones D., Karakas A. I., 
Koppen J., Tyndall A. A., Mohamed S. S., Rodriguez-Gil 
P., Santander-Garcia M., 2013, MNRAS, 436, 3068 
Miyaji S., Nomoto K., Yokoi K., Sugimoto D., 1980, Publi¬ 
cations of the Astronomical Society of Japan, 32, 303 
Mowlavi N., 1999a, A&A, 344, 617 
Mowlavi N., 1999b, A&A, 350, 73 
Norris J. E., 2004, ApJ, 612, L25 
Norris J. E., Da Costa G. S., 1995, ApJ, 447, 680 
Paczyhski B., 1970, AcA, 20, 47 
Paczyhski B., 1974, ApJ, 192, 483 

Pasquini L., Manas P., Kaufl H. U., Cacciari C., 2011, A&A, 
531, 35 

Piersanti L., Cristallo S., Straniero O., 2013, ApJ, 774, 98 
Piotto G., Villanova S., Bedin L. R., Gratton R., Cassisi 

S. , Momany Y., Recio Blanco A., Lucatello S., Anderson 
J., King I. R., Pietrinferni A., Carraro G., 2005, ApJ, 621, 
777 

Renzini A., Voli M., 1981, A&A, 94, 175 
Roeser M., 1975, A&A, 45, 335 

Romano D., Matteucci F., Tosi M., Pancino E., Bellazzini 

M. , Ferraro F. R., Limongi M., Sollima A., 2007, MNRAS, 
376, 405 

Sackmann I. J., 1980a, ApJ, 241, L37 
Sackmann I. J., 1980b, ApJ, 235, 554 
Shingles L. J., Karakas A. L, 2013, MNRAS, 431, 2861 
Shingles L. J., Karakas A. I., Hirschi R., Fishlock C. K., 
Yong D., Da Costa G. S., Marino A. F., 2014, ApJ, 795, 
34 

Siess L., 2006, A&A, 448, 717 
Siess L., 2007, A&A, 476, 893 
Siess L., 2010, A&A, 512, 10 

Siess L., Goriely S., Danger N., 2004, A&A, 415, 1089 
Stancliffe R. J., Jeffery C. S., 2007, MNRAS, 375, 1280 
Straniero O., Cristallo S., Piersanti L., 2014, ApJ, 785, 77 
Straniero O., Dominguez L, Cristallo S., Gallino R., 2003, 
PASA, 20, 389 

Straniero O., Gallino R., Busso M., Chieffi A., Raiteri C. M., 
Limongi M., Salaris M., 1995, ApJ, 440, L85 
Straniero O., Gallino R., Cristallo S., 2006, Nuclear Physics 
A, 777, 311 

Sweigart A. V., 1999, New Views of the Magellanic Clouds, 
190, 370 

Travaglio C., Galli D., Gallino R., Busso M., Ferrini F., 
Straniero O., 1999, ApJ, 521, 691 
Umeda H., Nomoto K., Yamaoka H., Wanajo S., 1999, ApJ, 
513, 861 

van Raai M. A., Lugaro M. A., Karakas A. I., Garcfa- 
Hernandez D. A., Yong D., 2012, A&A, 540, A44 
Vassiliadis E., Wood P. R., 1993, ApJ, 413, 641 
Ventura P., Carini R., D’Antona F., 2011, MNRAS, 415, 
3865 


Ventura P., D’Antona F., 2008, A&A, 479, 805 
Ventura P., D’Antona F., 2009, A&A, 499, 835 


APPENDIX A: EXAMPLES OF THE ONLINE 
TABLES 

In Table [XT] we show the first few rows of an example online 
table that contains the stellar yields of each chemical element. 

This paper has been typeset from a T[^/ DT|gX file prepared 
by the author. 
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Table Al. The first few rows of an example stellar yield table. Each table begins 
with a header describing the stellar model. 


# Mini = 3.00 Msun, Z = 0.0006, Y = 0.24, Mpmz = l.OE-03 Msun 
^ Initial [Fe/H] = -1.41, [alpha/Fe] = 0.0 

# Mremnant = 0.801 Msun, Myieid(all) = 2.199 Msun 


El 

z 

log e(X) 

[X/H] 

[X/Fe] 

Myield 

^yield 

P 

1 

12.000000 

0.000000 

1.361876 

1.441618E-I-00 

6.556237E-01 

he 

2 

11.097055 

0.167055 

1.528930 

7.210486E-01 

3.279209E-01 

c 

6 

9.231870 

0.761870 

2.123746 

2.961464E-02 

1.346825E-02 

n 

7 

7.679011 

-0.190989 

1.170886 

9.638053E-04 

4.383226E-04 

o 

8 

7.569763 

-1.160237 

0.201638 

8.565194E-04 

3.895308E-04 

f 

9 

5.327235 

0.907235 

2.269111 

5.818873E-06 

2.646327E-06 

ne 

10 

8.093986 

0.123986 

1.485862 

3.918737E-03 

1.782176E-03 

na 

11 

6.140176 

-0.099824 

1.262052 

4.578822E-05 

2.082372E-05 

# 

[Rb/Zr] 

[Is/Fe] 

[hs/Fe] [hs/ls] [Pb/hs] [C+N+Q/Fe] 

log e(CNO) 


-0.2027 

1.4595 

1.9514 0, 

,4919 0. 

9292 1.9422 

9.5376 



